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Abstract — In this paper the problem of assessing bounds 
on the accuracy of pilot-based estimation of a bandlimited 
frequency selective communication channel is tackled. 
Mean square error is taken as a figure of merit in channel 
estimation and a tapped-delay line model is adopted to 
represent a continuous time channel via a finite number of 
unknown parameters. This allows to derive some proper- 
ties of optimal waveforms for channel sounding and closed 
form Cramer-Rao bounds. 

Index Terms — Estimation, Fading Channels. 

I. Introduction 

Channel estimation plays a critical role in modern dig- 
ital communication systems, where receivers often need 
to acquire the channel state for each transmitted data 
packet. To facilitate channel estimation, pilot signals, i.e. 
waveforms known at the receiver, are usually embedded 
in the transmitted data signal [ 1 ]. In any application, it is 
important to devise pilot signals in a way that, for a given 
figure of merit, optimality or near optimality is ensured 
in a wide range of channel conditions. Important exam- 
ples of such a figure are represented by the Cramer-Rao 
bound (CRB) and the Bayesian CRB (BCRB), which 
limit the mean square error (MSE) performance achiev- 
able by any channel estimation algorithm. These bounds 
have been evaluated for a pilot-aided transmission in 
single-input multiple-output (SIMO) and multiple -input 
multiple-output (MIMO) block frequency selective fad- 
ing scenarios in J2j, under the assumptions that: a) 
the pilot signal is generated by a digital modulator fed 
by a sequence of pilot data; b) a symbol-spaced discrete- 
time model can be adopted for data transmission and, in 
particular, for the representation of a multipath fading 
channel; c) the tap gains of the channel model are in- 
dependent and identically distributed complex Gaussian 
random variables (this assumption is made in |f3] only). 

In this correspondence we revisit the problem of 
assessing performance limits on pilot-aided channel esti- 
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mation over a frequency selective channel, taking a novel 
perspective. In fact, we adopt a continuous time (instead 
of a discrete time) model for the overall description 
of a channel sounding system and adopt the MSE of 
the estimated continuous time channel impulse response 
(CIR) as a figure of merit. Then, we show that bounds 
for this figure of merit can be derived exploiting CRB's 
referring to the estimation of the tap gains of a tapped 
delay line (TDL) model of the communications channel. 
This sheds new light on both the achieveable limits 
and the properties of optimal waveforms for channel 
sounding; in particular, the role played by the properties 
of a continuous time communication channel in limiting 
the MSE performance in channel estimation is unveiled. 
This Correspondence is organized as follows. In Sec- 
tion |II] the model of a system for pilot-based channel 
estimation is described in detail and two figures of merit 
for channel estimation are defined. Two bounds on such 
figures are derived in Section [En] and are evaluated in 
Section [TV] for two different scenarios. Finally, Section 
IVl offers some conclusions. 

II. Signal and System Models 

In the following we consider the channel sounding 
system illustrated in Fig. [Q In this system, the transmitter 
sends a bandlimited real low -pass signal x(t) (dubbed 
pilot signal in the following), having bandwidth B and 
known to the receiver, over a frequency selective commu- 
nication channel characterized by its impulse response 
h(t) (or, equivalently, by its frequency response H(f)). 
Let r (t) = x (t)<g>h,B (t)+n (t) denote the noisy channel 
response to x(t), where ® denotes the convolution 
operator, n(t) is a complex circularly symmetric additive 
white Gaussian noise (AWGN) characterized by a two- 
sided power spectral density 2Nq and 



h B (t) 



B 



H(f)exp(j27rft)df 



(1) 



is a bandlimited version of h(t); note that hs (t) 
fully describes the noiseless channel behavior in the 
time domain for any input signal whose bandwidth 
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Figure 1. Channel sounding system: baseband model. 
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does not exceed B. The noisy signal r (t) feeds a re- 
ceiver which accomplishes ideal low-pass filtering (with 
bandwidth B), followed by sampling at a frequency 
f s = 1/T a = 2B, where T s denotes the sampling period 
(in Fig. [T]t n = nT s represents the nth sampling instant). 
We assume that the impulse response of the low-pass 
filter is g (t) = 2Bsmc(2Bt), so that its frequency 
response takes on a unitary value in the frequency 
interval (—B,B); then, the filter response y(t) is given 
by 



y(t) = x (t) ® h B (t) + w(t), 



(2) 



where w (t) is complex bandlimited Gaussian pro- 
cess having zero mean and a two-sided power spec- 
tral density S w (f) = 2Nq for |/| < B and zero 
elsewhere; note that its autocorrelation function is 
Rw(t) = 4NqB smc(2Br) and its average statistical 
power is a\ = R w (0) = 4NqB. Sampling y(t) generates 
the sequence {y n = y(t n )}, which feeds a channel 
estimator. This processes a finite subset of elements of 
{y n } to generate an estimate h B (t) of Kb (t). It is 
important to point out that: 

1. Any channel estimation algorithm assumes a spe- 
cific parametric representation of the communication 
channel. In the following, we adopt the well known 
tapped delay line (TDL) model for a bandlimited com- 
munication channel (H and assume & finite memory (i.e., 
a finite number of active taps); for this reason, h B (t) is 
expressed as 



M*) 



2B JT h B ,ismc\2B(t-^y\, 



(3) 



l=-Li 



where 



1 
2B 



Kb,i = wn^B (U 



lB hB \2B 



l -B! B B Hij)e A m w) di 



2B 



for any I and L\, L2 > (the overall number of active 



taps3 is L = L x + L 2 + 1). 

2. For a given sounding waveform x (t), a measure of 
the accuracy of the channel estimate hs (t) is provided 
by the MSE, defined as 

<"+oo 
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-B 



\eB,L(t)\ dt 



\E B , L (f)\ 2 df 



(4) 



where e B ,L (t) = h B (t) - h B (t) and E B:L (/) 
= H B (/) - H B (/), if the CIR h B (t) is modelled as a 
deterministic unknown function, and as 



£B,L 



1 
2B 
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2B' 



E 
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w,hi 
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\e B ,L(t)\ dt 



-B 



\E B , L (t)\ 2 df 



(5) 



if h B (t) is modelled as an unknown random process. 
Here, H B (/) (H B (/)) denotes the Fourier continuous 
transform of h B (t) (h B (t)) and Ex {•} denotes a sta- 
tistical average with respect to the random parameter X. 
Substituting ©in © yields 



Li 



yit) 



J2 h B ,ix(t- — j +w(t), 

=-Li 



so that the sample y n can be expressed as y n = 
y{nT s ) = y (^j) = Ya=-u h B,ix n -i + w n , where 
x n = x(t n ) and w n = w(t n ). In our system model, 
the channel estimator processes the set of N consecu- 
tive noisy samples {y n , n = 1, 2, ..., N}, i.e. the noisy 
vector y = [y\, y2, ..., 2/at] t , to generate an estimate 
hs = [h B - Ll , hB,i-L x , -, h B: L 2 ] T of the L di- 
mensional channel parameter vector h# = [hs,-^, 
^s,i-Li> ••■. h B ^L 2 ] T . This results in the estimated CIR 

'Note that the values of the parameters L\ and L2 (and, conse- 
quently, the value of L) should be large enough to ensure a good 
accuracy in the representation of the bandlimited CIR hs (t) and, in 
particular, to capture most of the energy of this signal. For this reason, 
such values mainly depend on the power delay profile (PDP) of the 
considered channel and are not necessarily equal (further details are 
provided in Section ITVt . 
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h B (t) = IB Zt- Ll " h B,i sine (2B (t - £)). It is easy 
to show that: a) y can be put in matrix form as 



y = Xh B + w, 



(6) 



where w = [wi, u>2, •••, Wjv] is a vector of inde- 
pendent and identically distributed complex Gaussian 
random variables (each having zero mean and variance 
a%, = ANqB) and X is a N x L matrix whose element 
on its i-th row and j-th column is JQj = Xi+j+Li-i 
(with i = 1, 2, ..., N and j = — Li, 1 — L\, ..., Li); 
b) thanks to the property of orthogonality of the sine (•) 
functions appearing in the channel model ©, the MSE 
© can be also expressed as 

B2 f o ^\ B2 



£B,L 



J2 E, 



B.l 



h 



B.l 



J2 M SE 



B.l 



(7) 

i.e. as a scaled sum of the MSE errors associated with the 
L channel taps (a similar expression can be developed 
for eb,l ©)■ In the following Section the problem of 
deriving bounds for the parameters eb,l © and e~b,l 
© is tackled. 

III. Evaluation of Performance Limits on 
Channel Estimation 

In estimating the vector h^ defined in previous Sec- 
tion, it can be modelled as a vector of unknown deter- 
ministic parameters or as a vector of random parameters 
with given statistical properties. In this Section we 
take into consideration both models, deriving some new 
bounds on the channel estimation accuracy. 

A. CRB-based performance limit 

In this Paragraph we focus on the class of un- 
biased estimators of the unknown deterministic vec- 
tor h.B and derive a lower bound for the param- 
eter eb,l ©■ To begin, we note that eb,l can 
be evaluated as eb,l = Y2i=-l var (^-B,z)> since 
^w{\h B ,i - h B ,i\ 2 } = veLT(h B ,i) (with I = -L u 1 - L ls 
..., L2), where var(X) denotes the variance of the 
random variable X. A lower bound to v&r(hB,i) for the 
above mentioned class of estimators is represented b 
the CRB O, which, in this case, can be expressed a 
vav(h B ,i) > [Jc 1 ( h B)] tl with I = -L 1 , \-L x , ..., L 2 , 
where 

2 The independence of noisy samples is due to the fact 
that E{wiw* k } = R w {ti - t k ) = W Bsinc(2B(ti - t k )) = 
4NoBsinc(l — k) = if I 7^ k. In other words, noise samples 
are uncorrelated and, being jointly Gaussian random variables, are 
statistically independent. 

3 Note that, to ease the reading, the indices of the rows and of the 
columns of Jc (hs) and J^ 1 (h^) range from —L\ to L2. 




=h B 
(8) 

with I, p = —L\, 1 — L\, ..., L2, is an L x L complex 
matrix, known as Fisher Information Matrix (FIM), 
/y(y;hg) is the joint probability density function of 
y © parameterized by the unknown (random) vector 



T 



h B and h B = \h B -u, b,B,\-L x , -, h B ,L^ 
(deterministic) trial vectoiQ. Then, the lower bound 



is a 



l=-Li 



1,1 



tr^- 1 (h B )) (9) 



can be formulated for £b,l, where tr(A) denotes the 
trace of a square matrix A. From the model © it can be 
easily inferred that, given he = he, y ~ CJ\f(/j,, R w ), 
where fi = Xh B and R w = cr^Ijy is the covariance 
matrix of w (Ijv is the N x N identity matrix), so that 
the element on Z-th row and p-th column of Jc(h#) 
can be expressed as (e.g. see @ Paragraph 2], Q rel. 
B.3.25]) 



[J c (h s )] /)P = 2Re 



dp, 



H 



8h 



B.l. 



+ tr r: 



-1 



5R V 
dh^.i 



R„ 



-Rv, 



dp 



9h B, P 

5Rw^| 



(10) 



with I, p = —L\, 1— L\, ..., L2, where Re(x) denotes the 
real part of a complex number x. It is easy to show that 
dp/dh Bp = (l/2)(l+j')[xi_ p , x 2 - p , ..., x N - p ] T , where 
Im(x) denotes the imaginary part of a complex number 
x. Then, substituting this result in (TTOl and keeping into 
account that dH w /dh^ = On (where Otv denotes the 
N x N null matrix) yields, after some manipulation, the 
expression 

1 f N 

[Jc(hfl)],, p = -2-Re \J2 

w lm=l 




(11) 



k=l—p 



4 The trial vector is used to indicate that the differentiation operation 
in the FIM definition is against a deterministic (versus random) 
complex variable. In particular, if f(fi) is some function of the 
deterministic complex vector /x = [. . . , Uj, . ..], then the usual 

applies. 



definition ^^ = ± ( dfM + i dfM 

2 I SRc{ Mi } + "'aim{ Mi } 



9/J-* 
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The last result shows that the FIM depends on the sample 
sequence {x^} of the channel sounding waveform x(t), 
but is not influenced by the parameters of the TDL 
channel model. We are interested in optimizing the lower 
bound (0 (i.e., in minimizing its right hand side) with 
respect to such a waveform. To tackle this optimization 
problem we assume that x(t) is a sample function of a 
bandlimited random process having the following prop- 
erties: a) it is wide sense stationary (WSS); b) it has zero 
mean and power spectral density (PSD) S x (f) > (= 0) 
for / € (-B,B) (/ £ (-B,B)); c) its autocorrelation 
function R x (t) tends to for r — > oo more quickly than 
1/t; d) it is ergodic in autocorrelation. These assump- 
tions entail that: 1) the sample sequence {x n = x(t n )} is 
a discrete-time WSS random process having zero mean, 
autocorrelation function R x [l] = R X (IT S ) and power 
spectral density 

S x (f) = Y, R x [k]e X p(-j2nfkT s ) 

k=— oo 
+oo 

= f s J2 sM-if.y, d2) 

l=— oo 

2) R x [l] decreases more quickly than \/l for I — > oo, so 
that the series X]£-oo l-^HI * s convergent; 3) {x n } is 
ergodic in autocorrelation. Under the above assumptions, 

the equality limjv-xx) jj J2k=i- r, x l x k+ P -i = R x[p ~ I] 
holds with unit probability (see (11 lb), so that for a finite 
(and large) N (i.e., when a large number of samples of 
the received signal is available for channel estimation) 
the element [3c (hfi)]; p (fTTl ) can be approximated as 



N 



[3c{h B )\ p = -^R x \p-l]. 



ot 



(13) 



The adoption of this approximation leads to a real sym- 
metric Toeplitz FIM; this implies that: a) any eigenvalue 
of J^ 1 (he) is always not smaller than inf (S x (f)) 
00 lemma 4.1], so that tr(J^ 1 (h B )) (see ©) grows 
unlimitedly as L — > oo (this means that, for a given N, 
as the number L of channel parameters to be estimated 
increases, the overall MSE diverges); b) the following 
asymptotic result holds [8. theorem 5.2c]: 



lim I t r(J c 1 (h B ))=T s ^ 

L— »oo J_j 



N 



/a/2 



df 



f s /2 s x (fy 



(14) 



since NS^f)/^ belongs to the Wiener class (i.e., the 
sum of the absolute values of the FIM diagonal elements 
remains bounded as L — >• oo; in other words (N/a^,) 

+oo 

Yl \Rxx[l]\ < oo), S x (f) is a real valued function 



and S x (f) > for any /. Then, from © and COD the 
lower bound 

.2 rS.ll df 

(15) 



lim £ £l L > ^ _ 

L^oo L ~ N J_ fs/2 S x (f) 

can be easily inferred. This result depends on the power 
spectrum S x (f), which can be optimized to improve 
the quality of channel estimation under the constraint 

J" Icy 

T s J_f 12 S x {f)df = P x on the average statistical power 
P x of {x n }. Applying the method of Lagrange mul- 
tipliers to this optimization problem leads to the con- 
clusion that the right hand side of ([T5l is maximised 
(under the given constraint) if S x (f) = P x for any 
/ G (—f s /2, f s /2), i.e. if the power spectrum of {x n } 
is uniform (equivalently, R x [l] = P x 5[l]); this occurs if 
(see (fT2l ) 



S x (f) 



u 





Pol 

f. 



% fe(- fs/2, fs/2) 
elsewhere 



(16) 
since x(t) is bandlimited to f s /2 = B Hz. It is important 
to note that, if the optimal power spectrum is selected 
for {x n } and the approximation ([T3l is used, ([T3l gives 
[3c (h B )] Ji? = (N P x /al)S[p-l] and the FIM 3 C (h B ) 
can be put in the form 



N P x 
Jc(h B ) = — 5— It, 



(17) 



so that 3 c l {h.B 
[3 c l (W 1 



{a 2 J{NP x ))l L and var(/», 



> 



N P* 



\ll - N P x - WWR' Where SNR = P *K 
is the signal-to-noise ratio, and the bound (© becomes 



£B,L > 



L 



Pb,l- 



(18) 



N -SNR 

This results evidences that, for a given SNR and a 
given number N of processed samples, an increase in 
the number L of significant CIR taps is expected to 
have a negative impact on the quality of CIR estimates. 
Finally, it's worth noting that the result expressed by 
(fTTl) is similar to that derived in [2 Paragraph 3.1] 
for channel estimation based on a training sequence 
that consists of a large number of uncorrelated channel 
symbols. In El Paragraph 3.1], however, a discrete-time 
communication model is assumed in the derivation of 
Cramer-Rao bounds. 

B. BCRB-based performance limits 

In this Paragraph we assume uncorrelated scatter- 
ing (US) and model the CIR hs(t) as a complex 
Gaussian process characterized by a zero mean (i.e., 
Rayleigh fading is assumed) and a PDP Ph(t) with 
I-^Ph{r)dT = 1. Then, we have that h B ~ CAf{0 L , 
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Rh), where Rh is the covariance matrix of h#; the 
element on Z-th row and p-th column of Rh is given 
by (see ©) 

E {h B ,l h% p ] = k{^J b b H(h)e^ l ^dh 

'B-f. 



dh) 



can easily be developed for Eb,l ©• To evaluate the 
right hand side of the last inequality, let us compute 
now the partial derivatives appearing in (122V It is easy 
to show that 



<9m/ hf 



1 



(2B) Z J h= - B 



f=-B-f 2 



J277 



C-p)/2 



(19) 

with I, p = —L\, 1 — L\, ..., L2, where 
R H (f) = E{H(f + f)H*(fo)} is the channel auto- 
correlation function (i.e., the inverse continuous Fourier 
transform of Ph{r)) and /o is an arbitrary frequency. 
Note that for I = p ([19]) yields 



9h B, P 

1 /a(hgR^h s 




Rh'hB 



(24) 



E 



{\h B ,l\ 2 } 



1 



(2B) 



13 



2 
-1/2 

j/=-l/2 



B-h 



R H (f)e^ l ^df 

f, --B Uf=-B-f 2 
1/2-2/ 

## (2Sx) e jM:c dx 

x=-l/2-j, 



Then, substituting d24J in d22J> yields 
3 h (h B ) = E hB J (R^hfl) (R-h 1 ^) 



a 



d/2 

r/y 

(20) 

Generally speaking, channel estimation algorithms can 
benefit from the availability of information about channel 
statistics to improve the quality of their CIR estimate. For 
such algorithms a lower bound to their MSE performance 
is provided by the BCRB \9, p. 957-958], which estab- 
lishes that MSE(h B ,i) > [3 B 1 {h B )] u with / = -L lt 
1 — L\, ..., L2, where 3b (he) is an L x L complex 
matrix, known as Bayesian Fisher Information Matrix 
(BFIM). The element on the Z-th row and p-th column 
of 3b (h) can be evaluated as |3, equ. 53] 

[3 B (h B )] l:P = [3 c (h B )] l:P + [3 h (h B \ p (21) 

where JcO^b) is the CRB FIM evaluated in the previous 
Paragraph and 

[•M h B)];, P - 

dln/ hs (he) fdlnf hB (h B 



(Rh 1 )" 



R 



h ' 



(25) 



since Rh is an Hermitian matrix. Like the CRB, the 
BCRB is influenced by the choice of the sounding 
waveform through Jc(h#) (see dUJ); in the following 
a uniform power spectrum is assumed for this waveform 
(see C[6]>). Then, substituting £T7]> and ([25]) in ([2T) yields 



J B (h B ) 



NF Xt ^ 1 

-T2^ IL + R h* 

u w 

( R" 1 

N • SNR I L + h 



iV-SNR 



(26) 



Unluckily, JbO^b) is not a Toepliz matrix and, as far as 
we know, no asymptotic result is available for the trace 
of its inverse. However, a simple expression for this trace 
can be derived if the Taylor series representation 



Jb 1 ^) 



00 



R,: 



N • SNR ^ 

fe=0 



iV-SNR 



(27) 



Ev 



dh* 
un B,l 



dh* 
B, P 




=h B 

(22) 



where f^ B (he) denotes the joint pdf of he- Like in 
the previous case (see © and ©) the bound 

Li 



can be adopted for 3 b x {\\b)\ this holds if the L eigen- 
values of the matrix (l/(iV • SNR))R hl 1 are distinct and 
their values are less than unityQ i.e. l/(iV-SNR- Aj) < 1 
(or, equivalently, Aj > N g NR > 0) for i = 1, 2, ..., L, 
where {Aj, i = 1, 2, ..., L} denote the (real) eigenvalues 
of Rh- In fact, this representation entails that 



£b,l = 2^ Ew,/ib,, ^ h,B,i - hs,i > 

I = — L/ 1 



> 



trjj^ 1 ^)} 



1 



iV-SNR 



X> 



R, 



k=0 



AT. SNR 



(23) 




J2 [3 B 1 (hB)} l>l =tr(3 B 1 (h B ))^h,L 



The eigenvalues of the covariance matrix Rh are always positive; 
this implies that the eigenvalues of the matrix R^ 1 are also positive. 
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Since Rh is an hermitian matrix, its inverse R h 1 can 
be factored as R^ = U S^U^ [10, p. 245, sec. 5.2], 
where U is a L x L unitary matrix (whose columns are 
the eigenvectors of Rh) and ~E = diag{Ai, A2, ..., Xl}- 
Exploiting this factorisation it can be easily shown that 



tr 



K ±L) > ■-»■■ 



-1 



N • SNR / 



N -SNR 

k L 



l-^—\T- (29) 



since tr{UDU fl } = tr{D} for any matrix D (this 
result is known as similarity invariance property of the 
trace operator). Then, substituting the last result in 
yields 



tr 



_. 00 L / 

k=0 i=\ v 



-1 



N • SNR • \i 



•SNR^l + ^O, 

y I 

A-/ AT . QMR _l_ J_ ' 



N 

L 



i=l ± ^ AT.SNR-Ai 



, N • SNR + T 

1=1 A 



(30) 



since we have assumed that l/(N ■ SNR • Aj) < 1 for 
i = 1, 2, ..., L. Finally, substituting (1301 in d23l yields 
the bound 



1 



^ " ^ iV-SNR + i 



fo.L. 



(31) 



i=i 



It is worth noting that this bound depends on the 
statistical properties of the channel through the eigen- 
values of the matrix Rh, whose structure is re- 
lated to the shape of R B (/) (or, equivalently, of 
Ph{r)). Let us try now to simplify this bound un- 
der the assumption that the bandwidth B of the 
sounding signal is substantially larger than the co- 
herence bandwidth B c of the communication channel 
{wideband channel sounding). In this case we have 
thag (see <T5)) lfsJ' B _ fa RH (/)exp (j2tt/^) df = 
Phils) ~ p h(0) for any f 2 G (-B,B), so 
that E(h B jh* Bk \ =i P h (0)/2B if I = k and 
= if I 7^ k. Then, the channel taps are un- 
corrected, R^ 1 = (2B/Ph(0))I L , and (see dH) 

6 This approximation is motivated by the fact that B c provides an 
indication of the width of Rs(f) (i.e., of the frequency interval over 
which Rn{f) takes on significant values). Then, if B > B c , the 
following integral is negligibly influenced by a change in the center 
(J2) of the integration interval. 
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B = l/r d3 


(1,5) 


(3,4) 


(1,6) 


(1,6) 


B = W/r ds 


(1,48) 


(33,33) 


(1,61) 


(1,63) 



Table I 

Values of the couple (Li, L 2 ) capturing at least 90% of 

THE OVERALL AVERAGE ENERGY OF h B (t). 



J B (h B ) = ( N ■ SNR + t^t) II, so that the bound 
(l23l) becomes 



£b,l > 



L 



± P { bI (32) 



N-SNR + 2B/P h (0) 

Note that 2B/Pf l {Q) >> 1 because of the assumption 
of wideband signalling over the communication channel. 
Therefore, a comparison of the last result with (ITFi ev- 
idences that, in this scenario, a significant improvement 
in the quality of channel estimates should be expected if 
the channel estimator is endowed with a knowledge of 
the channel statistics. 

Finally, we note that the result (|32l is substantially 
different from the BCRB evaluated in Appendix 
A], which refers to a discrete-time channel model in 
which the channel taps are independent and identically 
distributed random variables with a given pdf. 

IV. Numerical Results 

The bounds expressed by ([Tol l and (l23l (with Js(h) 
given by d26l )) have been evaluated for an exponential 
(E), a Gaussian (G), a uniform (U) and a truncated ex- 
ponential (TE) PDP [11], so that P h {r) = e ~ T ^ ds u(r), 



Ph(r) 
Phir) 



'/(2t 



u(t)-u(t-t m ) 
r (l-e- T A// T o)' 



Ph(r) 



u(T)-u(T-T d3 Vl2) 



r ds V12 

-t/to respectively, where u(r) 
is the unitary step function, r^ s is the rms channel 
delay spread, tm is the maximum delay in the TE 
PDP and tq is another time parameter depending on 
Tds (see [11, eq. (16)]). In our simulations the channel 
bandwidths B = 1/r^ and B = lO/r^ (wideband 
channel sounding) have been taken into consideration. 
In both cases and for each of the above mentioned 
PDP's we have evaluated the smallest values of the 
parameters L\ and L2 ensuring that the overall average 
energy Ya=~l x K {\hB,i\ 2 } (where E{\h B ,i\ 2 } is given 
by dlQl)) associated with the RHS of © is at least 
90% of the overall average energy of Kb (t) (see Table 
IJ). Then, on the basis of such values, the couples 
(Li,Li2) = (3,6) and (Li,!^) = (33,63) have been 
selected for B = I/t^s and B = W/rds, respectively, 
since they encompass all the cases of Table [U Fig. [2] 
(Fig. |3]) illustrates the bounds j3 B ,L CLB and J3 B ,L dHJ 
versus the SNR for B = l/r ds (B = W/r ds ) and all the 
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Figure 2. Performance bounds /3b, l GD an d Pb.l d23b versus the Figure 3. Performance bounds /3b,l J18t and /3b,l J23t versus the 
SNR for different PDP's in the case B - l/r ds . SNR for different PDP's in the case B = 10/r ds . 



considered PDP's. These results show that: a) indepen- 
dently of the bandwidth adopted for data transmission, 
the impact of the availability of a priori information on 
the estimation accuracy of a communication channel is 
significant mainly at low SNR's (where the terms {1/Aj} 
(PH , not included in ( fT8l ), yield a performance floor); b) 
the BCRB is negligibly influenced by the PDP type; c) 
there is a significant performance gap between the case 
B = 10/r^s and B = l/r^ (this is due to the fact that 
the overall number of channel taps to be estimated in the 
latter case is substantially smaller than that of the former 
one). Our simulations have also evidenced that: 1) in the 
considered scenarios an accurate approximation of (123) 
is provided by eq. (OH for both values of B; 2) eq. (l32l 
represents a loose bound for the case B = 10/rds- 

V. Conclusions 

The problem of assessing performance limits on pilot- 
aided channel estimation of a time-continuous frequency 
selective channel has been investigated. Novel bounds 
based on the CRB and the BCRB for TDL channel 
models have been derived and have been assessed for two 
different scenarios. The derived results shed new light on 
the achievable limits of pilot-aided channel estimation 
and the properties of optimal waveforms for channel 
sounding. 
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